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Abstract 

Weyl's expansion for the asymptotic mode density of billiards consists of the area, length, curvature 
On I and corner terms. The area term has been associated with the so-called zero-length orbits. Here closed 

nonperiodic paths corresponding to the length and corner terms are constructed. 

" . PACS numbers: 05.45., 03.65.Sq 

^ ! 1 Introduction 

f^ , The asymptotic mode density appears in various branches in physics. In a very early analysis of the density 

On ' of eigenmodes for a cavity with reflecting walls, Weyl proved that the leading term is proportional solely to 

ON \ the volume of the domain, and independent of the shape [y . Since then terms which account for the surface, 

(~| ' shape and connectivity of the domain were found to refine asymptotic expansions |g, |[ ^. For a billiard 

CIhI inside a simply connected domain of the Euclidean plane, the two-dimensional version of Weyl's expansion 

reads 
is ' AC 
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^ ' where Dirichlet boundary conditions are used, units are set at 2m = fi = 1, the boundary is characterized 

by the arc length s, and A, C, c{s), a^ are the area, the total length, curvature and inner angle of corner, 
k> " respectively. 

'V^ . Most studies of eigenvalue densities are based on the Green's function technique. Mathematical asymp- 

C^ I totic theory involves the Tauberian theorems from the theory of Laplace transforms. Balian and Bloch worked 

with energy Green's functions. Without requiring Tauberian theorems, they obtained a multiple reflection 
expansion for Green functions. The use of curvilinear coordinates was also included in their discussions. 

The area term has the simple meaning that the probability for a system to be in a particular subset of 
phase space is proportional to the volume of the subset. This term has been associated with the so-called 
zero- length orbits in the semiclassical theory g. 

The periodic orbit theory of Gutzwiller, known as the trace formulas, relates the fine details of oscillating 
density of states with classical periodic orbits g 0] . A trace formula has been derived for integrable systems 
by Berry and Tabor [gj. Other extensions of the Gutzwiller theory have been developed g, |l^. However, 
not much attention has been drawn to the role of nonperiodic orbits. Here we shall relate some nonperiodic 
orbits to Weyl's expansions. In Sec. II a family of closed nonperiodic orbits of single reflection is given to 
account for the length term. In Sec. Ill orbits contributing to the corner term are examined. Finally, some 
remarks are made in the last section. 

2 The length term 

Using the stationary method, Gutzwiller derived from Feynman's path integral the approximate propagator 
from r to r' in the time difference t Q 

Ksc{v, r'; i) = ^ (27ri)~^/V|detC|exp [iR{v, r'; t) - iM7r/2] , (2) 

cl.tr. 



where R is Hamilton's principal function, M the phase index obtained by counting the number of conjugate 
points along each classical path over which the summation is taken, N the dimensionality, and C the matrix 
of the negative second variations of R 

By converting time to energy, the propagator is transformed to Green's function as 

POG 

G{r,r';E) =i-^ dtK{r,r';t)exp{iEt). (4) 

Jo 

By using the propagator (0), after integrating over time again by the method of stationary phase, Gutzwiller's 
semiclassical approximation to Green's function is given by 

G.e(r, r'; E) - ^^^^^^^ Y. V^^xp [z^(r, r'; E) - zm^/2] , (5) 

^ ' cl.tr. 

where S is the action integral, n is again the phase index, i.e. the number of conjugate points, but obtained 
by varying the trajectory at constant energy E instead of time t, and 



i^ = (-ir^detC=(-ir-det^^J, (6) 

where in the last equation r^ is in the subspace transverse to the trajectory in the local coordinate system, 
and p_L is its conjugate momentum. Once Greeen's function is known, the density of states can be calculated 

by 

p{E) = --Im I drG{r,r;E). (7) 

When the semiclassical approximation is made, the density of states is expressed in terms of classical 
closed orbits. The zero-length orbits contribute to the average density of states. When the method of sta- 
tionary phase is again used for the integration over r, only the periodic orbits still remain in the summation. 

For a billiard, there is a family of closed orbits involving single reflection. Consider the simplest case when 
the boundary consists of the a;-axis, and the inside region is the upper plane. A closed orbit of the family is 
an orbit going from (x, y) (with y positive) to [x, 0) and then returning back to {x,y). The calculation of the 
Jacobian dr±/dp'j^ in the expression (H) for D is given in the appendix. From Eqs. (^), and (|4^), noticing 
that 

E = k^, L = 2y, q = 2k = L/t, d^S/dE^ = -L/{Ak^), 

we find 

D^l/{8ky), |dctC| = l/(4i2). (8) 

It is easy to verify that Hamilton's principal function R and the action integral S are 

R = yVt, S = 2ky. (9) 

From the approximation (0) for the propagator, we have 

i.eo = -^exp(4), (10) 

where subscript 'co' stands for 'closed orbits'. By using the formula for Hankel's function 



Hl'^\xz) = z" / exp 
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from Eq. m) Green's function is obtained as 
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GUr,r;E) = --H'o>{2ky). (12) 

By means of the relation (R) , its contribution to the density of states may be estimated by using 

l + /i + i^\ f 1 + II — v^ 



to be 
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which gives exactly the length term in Weyl's expansion m). 

If we estimate Green's function directly from Eq. (0), the result becomes 



G^„ = --4^exp(2ifcy), (15) 



which is consistent with the above Geo in the asymptotic approximation. However, its contribution to the 
density of states is C/{4V2Ett), which is different from the expected result by a factor of V2. This situation 
is similar to that of zero-length orbits. Since integration over small y contributes significantly, we have to 
use the uniform approximation with Hankel's function. 

When a bounce takes place at a position s of the boundary with curvature c{s), the closed orbit along 
the normal to the boundary still exists. Set the coordinate system with the origin at the bounce point and 
the j/-axis along the normal directing towards the inside of the bilhard. By using Eqs. (39) and (pO|), the 
counterparts of quantities (@) are found to be 



D = l/[8ky{l-cy)], \dctC\ = l/[4t^{l ^ cy)], (16) 

which gives the same length term. 

3 The corner term 

In the previous section we have considered closed orbits of single reflection. For a corner with an acute angle 
there is a family of closed orbits with double reflection. We shall first examine an acute corner, and then 
extend the analysis to an obtuse corner. 

3.1 Acute corners 

For a corner with an acute angle a we may construct a closed orbit of double reflection as follows. Let us 
denote by O the vertex of the corner, and by OA and OB its two sides. Suppose that the mirror image 
of OA with respect to OB is OA' , and the image of OB with respect to OA' is OB' . In this way we get 
the first and second images of the original corner, and a given point Q is then mapped to Q2 in the corner 
A'OB' , as shown in Fig. 1. The straight line QQ2, after mapping its segments back into the original corner, 
gives the closed orbit looked for. 

Denote by r = \OQ\ the length of OQ. The length of the closed orbit is then 

IQQ2I =2rsina. (17) 

From Eqs. (H) and (||), we have 

^, 1 / (rsina)^\ ,,„, 

i^co. = ^exp .^-^ , (18) 



which, from Eqs. (m and (M), leads to 



—-^S{E), (19) 
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where we have used formula (03) and have written E a.s E + ie to give 



' -V^ + ^-SiE). 
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The corner term in Weyl's expansion for the rectangular corner is estimated to be 1/16, in agreement with 
the value given by expression (|lS|). 

Expression ( |l9| ) is not the same as the corner term in Weyl's expansion (nh. A corner gives a correction 
to the length term due to the restriction on the domain of integration. Set the vertex O at the origin, and 
side OA along the x-axis. The change in the domain of integration 

/•oo /• /"y^ 

dx dy — > dx dy, with 7 = tana 
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gives the correction to the density of states 

Ep = — Im \ dy\ dx -H^o\^ky) = -—5{E), (20) 
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where the factor 2 in front accounts for the two sides, and use of Green's function (O) for single reflection 
orbits has been made. Thus, by combining the above two corrections, the total is 

Pco2 = ^f-^ + 2cota)<5(£;). (21) 

Stt \ sm a J 

This is not of the same form as that in Weyl's expansion. For a = (1/2 — 5)^ a little smaller than the 
rectangle, from Eq. (pTl) the lowest order is 1/16 + 5/8, while that from Weyl's expansion is 1/16 + 5(5/24. 
For very small a, compared with 7r/(24a) from expression (1^), the semiclassical value is 3/(87ra). The ratio 
is 9/7r2 w 0.912. 

In the above discussion we have considered only orbits of double reflection which hit the side OB first. 
An extra factor 2 should appear for those orbits to count the two different ways of selecting orbits according 
to which side is hit first. 

3.2 Obtuse corners 

For an obtuse corner the above closed orbits of double reflection do not exist. A natural way of continuation 
has to be found. For this purpose we use the folding property of the propagator 

A-(r.r';,-0 = /*"A-(r.r";*-OA-(r"..';<"-0. i^ 

to include two-piece closed orbits. We first examine the case of an acute corner. Let us consider a two- 
segment broken line from Q to Q2 via some mediate point Q' shown in Fig. 1, and set t — t" ^ t" — t' — t. 
Suppose that the polar coordinates of these three points are (r, ^1), (r2,02) = (r, 2a + 9i) and (ro,0o), 
respectively. In similarity to Eq. (|lO[), we may derive the propagator for each segment. By means of the 
folding property, we find for the broken path 



K{Q, Q2\Q') = ^-^ I dOo exp 1^ [2r^ + 2rl - 2rro cos(0o - ^1) - 2rro cos(0o - ^2)] | 



(23) 



where the integration domain for ^o is determined by the constraints < 9q < 3a, < 6*1 < a, |6'o — 6'i| < tt 
and 16*0 — 6*2] < tt. At the limit of the rectangular corner, propagator (p3|), after integrating over 6*1, reduces 
to 



j deiKiQ,Q2\Q') = ^KiQ,Q2\Q')^ 
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exp U- . (24) 



Compared with the expected form (p^), K{Q, Q2IQ') equals the half of Kco2- It can be verified that if we 
use the method of stationary phase approximation in the Cartesian coordinate system for the integration 
involving Q' the propagator (p3) would revert to (O) exactly. However, if we count the two ways of selecting 
orbits, the result here turns out to be better. 

Generally, we may approximate the propagator between a reference point and some mediate one, say 
Q and Q' in Fig. 1, by the semiclassical propagator involving classical paths with to 2 bounces. Besides 
images Qi and Q2, point Q has two more images Q-i and Q-2 obtained clockwisely, as shown in Fig. 1. 
The paths which contribute to the semiclassical propagator from Q' to Q arc straight segments Q'Q, Q'Qi, 
Q'Q2, Q'Q-i and Q'Q-2- The last two correspond to paths hitting side OA first. Similarly, by considering 
images of Q', paths contributing to the propagator from Q to Q' can be found. By means of the folding 
property, combination of these two propagators gives an approximate propagator from Q to Q itself, which 
includes K{Q,Q2\Q') as a part. 

Calculation in the Cartesian coordinate system for the rectangular corner is rather easy. In this case 
both Q2 and Q-2 are in the third quadrant, while Qi and Q-i are, respectively, in the second and fourth 
quadrants. Each path arises as a length square in the exponent of the expression for propagator. Denote 
by {x,y) and (a;o,2/o) the coordinates of Q and Q' , respectively. There is a correspondence between closed 
paths and length square sums [{x ± xq)'^ + (y ± j/o)^] + [{x ± a^o)^ + (2/ ± Vo)^]- Here the first square bracket 
corresponds to the path from Q to Q' , and the second to the return path. A plus sign between x and xq 
indicates a bounce on side OB, while that between y and yo indicates a bounce on side OA. Thus, each 
path has its 'four signs' signature, which is the four signs appearing in the square sum. For example, the 

closed path without any bounces may be marked as , which contributes to the area term. It can be 

verified that each of the paths — I h and H 1 — results in 
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for the level density, while each of the paths H , — \ , \ — and h gives 

1 £ 1 
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Closed paths with a total of 2 bounces, besides paths — \ h and H \ — , arc H h, — hH — , +H 

and 1-+, each of which contributes 5(£')/64. Paths with 3 bounces are + + H — , + H h, H h+ and 

— 1-++, each of which contributes —S{E) /32Tr. The only path with 4 bounces is + + ++, whose contribution 
is S{E)/{16tt'^). The total contribution of these 16 terms recover exactly the area and length terms, and 
gives the corner term as 

We see that the paths with a single bounce on both sides gives the main contribution to the corner term. 

Thus, for an obtuse corner, although closed classical orbits of double reflection generally do not exist, we 
may still calculate the corner term from the two-piece paths of just 2 single reflections on both sides, which 
are made either by only one piece, or by each of the two pieces. Due to the cancellation among terms for 
the folding propagator the approximation keeps the main contribution. However, now the integrals involved 
cannot be estimated analytically, and numerical methods have to be used. 

4 Discussions 

In the above we have examined the role played by closed classical orbits with single and double reflection on 
the boundary for billiards. Here we make some concluding remarks. 

1. So far, we have considered only the length term and the corner term, leaving the curvature term 
untouched. The same sign for both Dirichlet and Neumann conditions gives us a hint that orbits of 
double reflection is dominant. In principle, we can consider the contribution from two-piece closed 



pathes, as we did for obtuse corners. Unfortunately, even the simple case of a circle does not admit 
simple expressions. One way to get round the intricate situation of a circle is to consider its inscribed 
polygons g. 

2. A way to derive the propagator for a corner is to use curvilinear coordinates. A corner can be 'flattened' 
by introducing the tranformation from {x,y) — > {u,v) defined by 

,_ , . ,_ N -12 2 2 2 7^ TT _ a , , 

X — r cosHip), y — r smHip) , with r ^u +jv , tamp— — , •y — —, •y——, (25) 

u a t: 

where a is the inner angle of the corner as before. The Jacobian of this transformation is equal to 
one, which is made to keep the area and length terms unchanged. The transformed Laplacian can be 
derived as 

A = dldUl' - 1) {dl - pi - pi - ^9„a. - J9„ - ^d?j . (26) 

We now regard the product with the factor (7^ — 1) as the perturbation to d'ldl = Aq. The factor 
(7^ — 1) is indeed tempting if one notes that the corner term is 

By means of the perturbation expansion for the propagator the corner term can be obtained. 

3. Only Dirichlet boundary conditions have been considered in the above. The extension to Neumann 
conditions is rather straightforward. Since there seems to be no general formula for the corner term 
at Neumann conditions in the literature [^, a semiclassical estimation can now be made using our 
approach. 

4. The closed paths considered above have a zero limiting length. There are other closed orbits of a 
non-zero limiting length. For example, there is a continuous family of closed orbits from the diameter 
orbit to the equilateral triangle orbit in a circular disk. Any member of the family is an isosceles 
triangle with one vertex inside the circle. The role played by such orbits is worth examining. 

5. We may extend our analysis to include connectivity and higher dimensionality. 
Some problems are under study. 



Appendix: Jacobian dr^/dp'^ 
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For a billiard inside a simply connected domain of the Euclidean plane, we may derive the Poincare map 
from bounce to bounce in Birkhoff coordinates (s, t;), where v is the component of the velocity in the tangent 
direction to the boundary right after reflection, and s the arc length along the boundary. Since the absolute 
value of velocity is conserved for a billiard, we may normalize the velocity as a unit vector and then let 
V £ [0, 1]. The linearized Poincare map from (si, vi) to (s2, V2) can be expressed as [^ ^ |l3[ 

M(12) = ( (^i2Ci-^i^)M± -h2/vi^V2^ \ (28) 

\ CiV2± + C2Vi^ - I12C1C2 (/12C2 - W2_l)/w1_L / 

l/t'2± \ / 1 \ / -1 -1,2 \f 1 \ / t-i^ \ 

V2± J \ -C2/V2± I J \ -I J \ -Cl/Vn_ I J \ 1/VH_ J ^ ''' 

where I12 is the length of the chord joining si and S2, v± is the normal component of the velocity and c the 
curvature of the billiard boundary. It is often useful to know the Jacobian matrix d{s,v)/d{£,,K), where £, 
and K are perturbations in displacement and velocity at a given point O on the straight line joining si and 
S2 along the direction perpendicular to the path. Without loss of generality we may choose the coordinate 



system with the origm at O, and the j/-axis along the path from si to S2, as shown in Fig. 2. (In the figure 
si and S2 are marked as s and s', respectively.) It is obvious that at si 

V^^ =0, Vy^ 1. (30) 

Denote by a dot the derivative with respect to the arc length s. The tangent and normal unit vectors are 
ti = (ii,2/i) and n — {~yi,xi), respectively. This implies that 

Vi ~ Xi and Wij^ ~ yi. 

Assume that a perturbation {Ss,Sv) from si to s = si + Ss results in the perturbation at O, (^, k). Up to 
the lowest order, we have 

X = x{s) — xi + xi5s = vij^Ss and y = y{s) ~ yi + yi5s w j/i. (31) 

It can be seen that Svix ^ Vx — vix ~ k- From Eq. ( pO[ ), the relation Vx{Svx) + Vy{5vy) = implies that 
6viy = 0. That is, at s we have Vx — n, and Vy — viy = 1. Using Vx/vy — {x — ^)/y, we find 

vi±6si-£,=yiK. (32) 

Similarly, the relation v — VxX + Vyi/ leads to 

Svi — vi±K + ySs — vi±K + ciVi±Ss, (33) 

where we have used the curvature formula y — ex. Equations (p2|) and (BS) can be written as 



(5si \ _ / l/ui± yi/vi± \ f i 
Svi J \ ci vi± + ciyi 1 \ K 



(34) 



The above transform matrix may be written as 

/vi± 
vi± J \ ci/vi± 1 / I 1 



A.=('/^ :)(J... ")(i ",')■ (35) 



T-i _ I 1 -2/1 \ I 1 \ I vij_ 



Its inverse 

''■ - ■'.v n i" J I -/.'.x 1 ) [ T 1/;. ) • p« 

describes the tranformation from {6si, Svi) to (f , k). 

Along similar lines we may derive the Jacobian matrices between {5s2,Sv2) and (^, k). There is a main 
diS^erence. When following the above derivation for S2, we use Vx and Vy of the velocity right before a bounce. 
By taking this into account, vj_ in the above formulas has to be replaced by —v±. For example, 

Ss2 \ _ f -l/v2i_ -y2/v2± \ f ( \ (37) 



SV2 J \ C2 -V2± + C2y2 / \ « 

and 

Noting that /12 == j/2 — J/i, we can verify that 

M{12)^J,^{s2)JUsi). 

For an orbit which starts and ends inside the billiard, respectively, at (ro,po) and (rt,pt), and makes 
successive bounces at si, S2, ..., s„ in between, we have 

^* )^m( ^° ), M^JUsn)M{n-l,n)---M{l2)Js^{si). (39) 



The Jacobian d£,t/dKo is determined as 

dCt/dKo = kdr^/dp'^^Mi2, (40) 

where k is the absolute value of the conserved momentum. When all the bounces happen at straight segments 
of the boundary with curvature c = 0, matrix M is significantly simplified. In this case, using expressions 
(ll) and (H), we have 



where L is the total length of the orbit. Thus, from Eq. (|40|) we have 

dr^/dp'^ = {-iTL/k (42) 

for the Jacobian. 
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Fig. 1 Images of a given point Q in a corner. Some mediate point Q' is used to construct folded paths. 
Fig. 2 Perturbation of a path. 
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